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Abstract 



J> ■ We re-examine the dynamical stability of the nakedly singular, static, spherically 

symmetric solutions of the Einstein-Klein Gordon system. We correct an earlier 

proof of the instability of these solutions, and demonstrate that there are solutions 

ON ' to the massive Klein-Gordon system that are perturbatively stable. 
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1 Introduction 



We are re-considering the dynamical stability of the non-trivial static, spher- 
ically symmetric solutions of a real, self-coupled, self-gravitating scalar field. 
The case where the scalar field is massless and the potential vanishes has ana- 
lytically known static solutions [13] which, as claimed by Jetzer and Scialom [4] , 
are unstable against linear perturbations of the system. This is in accordance 
with the cosmic censorship hypothesis [12] and the results of Christodoulou's 
study of this system (beginning with [2]) which led to the conclusion that 
the set of initial data that collapse to these nakedly singular solutions is of 
measure zero [3]. 

That the massive scalar field has quite a different flavour is indicated by the 
presence of periodic in time, soliton solutions [8] and the richer structure re- 
cently discovered in collapsing solutions that are near the threshold of black 
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hole formation [1]. As one further piece of the puzzle, we show here that there 
are nakedly singular solutions to the massive scalar field equations which are 
stable against spherically symmetric perturbations. This is in contradiction 
to the results reported by Jetzer and Scialom [4]; however, their calculations 
contain a flaw that leads to the erroneous conclusion that all of these nakedly 
singular solutions are perturbatively unstable. The results from the corrected 
version of their variational bound are consistent with the explicit determina- 
tion of the perturbative mode. 

The system of equations that we are considering is straightforward to derive. 
We will state here the basic equations and conventions we use. We consider 
spherically symmetric spacetimes with the metric in the diagonal form [y = 
v(t,r) and A = X(t,r)) 

g = diag(e", -e\ -r 2 , -r 2 sm 2 (9)), (1.1) 

dynamically coupled to a scalar field with Lagrangian 

where the potential and its variation are given by 

V[<t>] = \m 2 <p 2 + \k<P\ StVty] = m 2 <f) + n<f. (1.2b) 

Note that we have chosen 167rG = c = 1 and <fi is taken to be dimensionless; 
this implies that m has dimension of an an inverse length and k an inverse 
length squared. 

Einstein's equations may be written = ^(d^[(p]d u [(f)]— g^ u V[(p]), from which 
we find 

Roi = $dt[\] = Idt^dM, (1.3a) 
e-^Roo + e" A i?n = \e~ x d r \u + A] = \e~ v {d t [<t>\f + | e - A (9 r [0]) 2 , (1.3b) 
R22 = ~e~ x d r [v — A] + 1 — e~ A = ±rV[0]. (1.3c) 

The variation of the matter action with respect to results in the wave equa- 
tion 

e-^ 2 [0]-|e"^[z/-A]9 t [0] 

- e- A ,9 2 [0] - e- A (f + \d r [v - \\)d r [<j>\ + = 0. (1.3d) 

In Section 2, we present the numerical solutions of the static limit of (1.3), 
with the Appendix devoted to the analytic solution of the static equations 
with vanishing potential. In Section 3 we set up the perturbative stability 
problem and present the numerically generated ground state energies of the 
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perturbative mode, indicating that the spacetimes with vanishing potential 
are unstable but that there is a branch of the spacetimes with non-vanishing 
scalar field mass which is perturbatively stable. In Section 4, we examine the 
variational method introduced by Jetzer and Scialom to put an upper bound 
on the ground state energy, showing that the corrected method is consistent 
with the results of Section 3. 



2 The Static Background Spacetimes 

Throughout, we will make use of the dimensionless radial coordinate scaled 
by the asymptotic Schwarzschild mass parameter M s of the system: x := 
r/(2M s ). In the case of a non-zero inverse length scale m or ^J~k~ we will define 
the dimensionless quantities m := 2M s m and k : = AM 2 k. We then re-write 
the static limit of (1.3) in a form appropriate for the numerical integration 
implemented below: 



where we have written z/ , Ao and 0o to represent the static background fields. 

We will employ the simple shooting method [9, Section 7.3.1] on 0o to numer- 
ically generate a solution to (2.1). Initial values for exp(z/ ), exp(Ao), 0o and 
<9x[0o] are given for small x and (2.1) are integrated outward using a fourth- 
order Runge-Kutta integration scheme [9, Section 7.2.1]. Assuming that both 
of exp(z/ ) and exp(Ao) vanish as some positive power of x as x — * 0, we find 



(ra+ := 2( V / TTF + 1 + k 2 )/k 2 E [2, oo] and x k := x/VTTW) 

e u ° ~ Fx n k + ~ 2 , e x ° ~ Gn\ x n k + , O ~ -k(n+ - 2) \n(x k ) + B, (2.2) 



independent of the values of m and k. The parameterization in (2.2) is moti- 
vated by the analytic rh = k = solution, which from (A. 8a) has F — G — 1 
and B = 0. (Note that when the potential vanishes the field equations are 
invariant under O ~~ * <Po + constant. This must be taken into account when 
applying the shooting method, however it is sufficient to check that B = 
and F — G — 1 reproduces the analytically known solution.) 

The small- x behaviour of these solutions indicates that there is a curvature 
singularity at x — 0, showing up in the Ricci scalar as R oc 1/x 4 . Further- 
more, by considering radial null geodesies emanating from x = (affinely 
parameterized by A and with tangent u), we find \ 2 R(u,u) ~ l/n^_; the finite 
limit of this as A — > indicates that there is a strong curvature singularity at 



d x [e Xo ] = e A °{f (cU0 o ]) 2 + ~ X {1 ~ e A °) + f e x °V[<p }}, 
d x [4*] = e^{| (dM) 2 ~ |(1 " e Ao ) - fe Ao \/[0 o ]}, 
d 2 M + (| + id x [u - Xo])d x [<j> ] - e^d^ViM = 0, 



(2.1a) 
(2.1b) 
(2.1c) 
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x = by an argument of Clarke and Krolak (see for instance [10]). For the 
solutions considered herein, it is straightforward to show that outgoing radial 
null geodesies are well-defined and reach x = oo in an infinite amount of afline 
time, as do outgoing radial timelike geodesies with sufficient energy. 

The large-x behaviour of 0o is determined by requiring that the metric com- 
ponents be of the asymptotic Schwarzschild forms 

e"° ~ 1 - i, e Ao ~ 1 + \. (2.3) 

Inserting these into (2.1c) we find (keeping only the asymptotically dominant 
terms and noting that we require that 0o —> as x — > oo) 

^[0o] + |4[0o]-m 2 0o~O. (2.4) 

From this, the asymptotic form of 0o is determined: 

\ A/x for rh = 

[Ae~ mx /x form^O 

and from (A. 5) and (A. 6) we know that for rh = k = we 

Given input parameters rh and k, the system (2.1) has a unique solution once 
the initial data k, F, G and B are chosen. Although the constant F may be 
freely chosen by rescaling the time variable (and is chosen by requiring that 
exp(z/ ) — > 1 as x — > oo), we must still determine fc, G and B. For each value 
of k that we consider, we avoid shooting on both of G and B by proceeding as 
follows: We choose values for rh and k and assign F — G — 1 initially, shooting 
on 0o to determine 5. This will result in a solution with asymptotic behaviour 
exp(Ao) ~ (1 + a/x) and exp(z/o) ~ /o(l — a/x). Re-scaling x by x — > ai and 
renormalizing exp(z/o) results in a solution with the correct asymptotic form 
for the metric functions, while the x ~ forms become 

e*° ~ /o-^ax^- 2 , e Ao ~ n\{ax k ) n + . (2.6) 

In addition, the constant B determined by shooting is shifted by B — > 5 — 
A:(n + — 2) ln(a). Most importantly, we have rh — > am and k — > a 2 /t, and we 
cannot treat rh and /c as fixed input parameters. Note that any solution that 
satisfies 0o ~^ as x — > may be rescaled to one that has F = G = 1, and 
therefore we are exploring the entire class of solutions in this manner. 

We have found that although k, F, G and B uniquely determine a solution 
of (2.1) and k and B parameterize the space of solutions with asymptotic form 
given by (2.3), for the k = solutions this does not uniquely determine rh. 
We show in Figure 1 two solutions with k — 3 and rh ~ 0.210 as well as the 
k = 3, rh = k = solution. From Figure 2(b) we see the following feature: 



(2.5) 
must find A — k. 
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above a minimum rh there is a stable and an unstable solution labeled by the 
same values of k and rh. (A similar result holds for the rh = solutions with 
k^O not considered in this work.) 
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Fig. 1. The background metric functions for three k = 3 solutions: the solid line 
corresponds to the unstable rh = k = solution, the dashed line the unstable k = 0, 
m 0.210 solution, and the dotted line the stable k = 0, rh 0.210 solution. 



3 Dynamical Stability 

We turn to the investigation of the dynamical stability of these spacetimes, i. e., 
whether they are stable against linear perturbations. Two different methods 
will be employed: a variational approach which casts the equation for the 
perturbative modes into the form of a stationary Schrodinger equation and 
uses a one-parameter family of wave functions to put a least upper bound on 
the lowest eigenvalue, and a shooting method to numerically generate both 
the lowest energy mode as well as the eigenvalue. 

In either case we make use of the single Fourier-mode expansions of the per- 
turbing fields &sv(t,x) = vq(x) + vi(x) cos(otf), X(t,x) = Xo(x) + X±(x) cos(u)t) 
and (f)(t,x) = (f>o(x) + (f>i(x) cos(ut)/x. The first-order perturbation equations 
from (1.3) for the metric functions give 

Ai = \d x [<f> Q }<f>i, d x [ui] = -|^[0 o ]0i - £3c[0o]<£i + i<U0o]dx[0i], (3.1a) 
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and the equation for 0i 



dl[<t>i] + \d x [u - \o\dM - e A °-^y o N0i + ^ 
where 
V [x] : 



o, 



{hd x [v Q - A ] - |(^[0o]) 2 - HWoD^h - A ] 
+ e A °(m 2 + 3^) + a;e Ao ^[0o](m 2 o + 



(3.1b) 



(3.1c) 



From the small-x behaviour of the background fields (2.2) we find 
d 2 x [<j )l ]-(l/x)d x [h} + (l/x 2 )h ~0, 



(3.2) 



from which we have that <pi ~ x(a + 61n(x)). We will restrict ourselves to 
b = so that the perturbation is finite at x — 0, however all of the solutions 
considered here appear unstable against singular perturbations with 6^0 
and therefore ln(x) behaviour at the centre of symmetry. Once the arbitrary 
amplitude a is chosen, (3.1b) is integrated numerically, varying the value of uj 2 
and requiring that <f>i — > as x — > oo. In Figure 2(a), we display the results of 
this procedure for the rh = k = solutions, giving the value of uj 1 for lowest 
energy mode. As we see, Co 2 < 0, indicating that these solutions are generically 
unstable. As discussed earlier, there are two solutions with given values of k 
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(a) rh = k = 0. 



0.05 0.1 0.15 0.2 0.25 



(b) k = for (from right to left) k 
3, k = 4, k = 5. 



Fig. 2. In both figures the solid line indicates the exactly determined eigenvalue 
from numerically solving (3.1b) and the dashed line the variational bound generated 
from (4.16). The diamonds indicate the solutions studied in Figures 1 and 3. 

and rh and, as shown in Figure 2(b), there is an unstable branch (uj 2 < 0) and 
stable branch (Co 2 > 0). 

In generating the perturbative ground state and eigenvalue for the rh ^ 
solutions, the presence of an exponentially diverging solution of (1.3) does not 
allow us to generate the background solution to arbitrarily large x. However 
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we shall see from (4.4) and (4.6) that for large enough p (and Co 2 < m 2 ) we will 
always have c^[0i] > 0, and therefore if 9 x [0i] > then the perturbation must 
diverge. When shooting on u 2 to determine the ground state, it is therefore 
sufficient to check whether (pi passes through zero (in which case u 2 is too 
large) or 4>\d x [(j)i\ > for large enough x that V — Co 2 > (in which case u 2 
is too small). 



4 The Variational method 



We now consider the variational method of generating an upper bound on 
u 2 , introduced by Jetzer and Scialom [4]. Although we find no fault with the 
method in principle, we show that they made a technical error in its appli- 
cation. Using their erroneous bound, one would find that these solutions are 
generically unstable, whereas the direct numerical determination of the ground 
state and eigenvalue in Section 3 contradicts their bound and furthermore in- 
dicates that there are solutions to the massive scalar field equations that are 
perturbatively stable (as shown in Figure 2(b)). 

Following [4] we introduce a radial coordinate p that eliminates the second 
term in (3.1b) 

p= f X dx e -| (iy °- Ao) , (4.1a) 
Jo 

chosen so that p = at x = 0. From (2.2) and (2.3) respectively we find that 



p/Fx 2 n + /{2VTT¥) as x - . ^ 
x + \n(x) + constant as x — > oo 

p is numerically determined by making use of the x — > result and integrating 
outward. 

This coordinate transforms (3.1b) to the stationary Schrodinger-like equation 

- d ^ + V[x{p)]<p 1 = ^ l . (4.2) 

We consider perturbations in the real Hilbert space L 2 (dp, (0, oo)), i.e., real 
functions on the half-line with the norm 

(0,0') = f°° dp<t>{p)<t>'{p) = r dxe-^-^mx^W)]. (4.3) 
Jo Jo 
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From the asymptotic and small-x forms given earlier, we find 



V[ P ] 



-V(V), v[ P ] 



and we therefore introduce 



1 



so that (4.2) is written as 



— 1/p 3 for rh = k = 
3kA 2 /p 2 for m = 0, k ^ , 
for m 7^ 



(4.4) 



m 2 



Ap< 



+ V[x(p)} : 



H <l> 1 + V[x(j>)]<f> 1 = u> 2 <l> 1 . 



(4.5) 



(4.6) 



From this form it is worthwhile to make a few observations. If we consider 
regions where V « Vq — constant, then for a) 2 > Vo we have the oscillating 

Bessel function solutions (defining Q := \J\Cj 2 — V \) 



0x « Ay/pJ (n P ) + B^pY (Clp), 



(4.7) 



and for cZ> 2 < Vo replace J and Yo with 7 and K to find the exponentially 
increasing and damped solutions respectively. From this we deduce that Co 2 < 
Hindoo V since otherwise from (4.7) <pi would not vanish as x — > oo, and 
from (4.4) we therefore find that the solutions with rh = can have no stable 
modes; we will not consider stationary perturbations (a) = 0) in this work. 
This has been verified for the rh = k = solutions as shown in Figure 2(a) as 
well as for the rh = and k ^ solutions. (Note that one must be careful in the 
latter case since for arbitrarily small and positive uj 2 , <pi will begin oscillating 
at arbitrarily large x even though it will appear to have exponential behaviour 
prior to that.) This does allow the possibility that the rh ^ solutions are 
stable, and from Figure 2(b) we have seen that there is a branch of solutions 
for which this is realized. 



From the work of Narnhofer [6] we know that the operator H Q is not self- 
adjoint on the intersecting domain £>(— <9 2 ) fl V(— l/(4p 2 )), with deficiency 
indices (1,1). A family of self- adjoint extensions is determined by extending 
this domain to include the solutions of 

-d 2 p [^ ±l ] - l/(4p 2 )^ = ±# ±i , (4.8a) 

which in this case are 

^ = ^p-H^ipe^ 4 ), ^ = ^ = y/pH^ipe-**'*), (4.8b) 

where and are the zeroth-order Hankel functions of the first and 
second kind respectively. 
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The self-adjoint extensions are parameterized by the real angle r, and extend 
the operator to act on linear combinations of the real functions 

^ T := \{e iT i>i + e-*^) = cos(r)^) - sm(r)3(^), (4.9) 

where 3? and 3 represent the real and imaginary parts respectively. (Note that 
we consider self-adjoint extensions on a real domain only since the scalar field 
is real.) We write the operator H t which acts like H on the extended domain 
£>(#o,r) = V(-d 2 p ) n £>(-l/(4p 2 )) + {^ T }, and using (4.8) acts on ^ r as 

#o,t*t = W T i>i ~ e~ iT iJ-i) = - sm(r)&(^) - cos(r)9f(^). (4.10) 
We will require the following integrals: 

/ dp K(^) 2 = / dp 3(>i) 2 =:a^ 0.1592, 

JO JO 

/ dp K(^)9f(^i) =:b « -0.1013, 

JO 

and using (4.9) and (4.10) we find 

(* T ,* T ) = (a-6sin(2r)), (* T , H , T ^ T ) = -b cos(2r). (4.12) 

This latter result corrects the error in the work of Jetzer and Scialom who have 
assumed that it vanishes (see the paragraph following Equation (28) in [4]). 

From the asymptotic forms of the Hankel functions, we find the behaviour of 
ty T for large and small p respectively to be 

^ r ~ y^re"^ sin(p/ V2 + tt/8 + r) , (4. 13a) 

^ r ~ v^cos(r)/2 - (2^/tt) sin(r)( 7E + ln(p/2)). (4.13b) 

At this point we note that only the r = extension (chosen implicitly in [4]) 
corresponds to a perturbative field that is finite at p = and, consistent with 
the choice made in Section 3, we restrict ourselves to this case in what remains. 

From the form of the differential operator Ho,o, f° r an Y t wo functions G 
V{Hq q) we define ij)p{p) '■= ip(f3p), and it is straightforward to prove the scaling 
property: Hopipp) = Hoflip). We therefore consider the collection of 
functions 

*oAp) ■= MPpX ( 4 - 14 ) 

as the variational family, varying (5 to get a least upper bound on the ground 
state energy. Using this scaling property and the results (4.11) and (4.12), we 
find 

(tf o,/j, #o,o*o,/3) = -bp, (*o,/3, *o,/3) = a/0, (4. 15) 



(4.11a) 
(4.11b) 
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and taking the expectation value of (4.6) leads to the bound 



Co 2 < /3((*o,/3, VVqj) - bf3)/a = : max(w 2 ). (4.16) 

The inequality (4.16) should be compared with Equation (30) of [4] which 
is lacking the —b(3 2 /a contribution. As we see from Figure 3(d), this term is 
necessary for the consistency of the variational result with the numerically 
generated eigenvalue, as it dominates for large (3. Naively we expect that very 




(a) m = k = 0-unstable. (b) k = 0, m = 0.210-unstable. 




(c) k = 0, rh = 0.210-stable. (d) 



Fig. 3. On the first three figures appear the effective potential (solid line) and the 
perturbative ground state (dotted line) for the solutions that appear in Figure 1. 
The final figure shows the behaviour of max(cl> 2 ) (4.16) with respect to (3 (labeled as 
in Figure 1): the curves with an obvious minimum correspond to the correct form 
of the bound (4.16), and those that diverge negatively are derived from (4.16) by 
dropping the second term and corresponds to that of Jetzer and Scialom. Note that 
this final form clearly contradicts the results shown in Figure 2. 

small values of (3 will result in a bound of Co 2 < m 2 since in these cases the bulk 
of the support of ^o,/3 is shifted to large p and the asymptotic behaviour of the 
V (4.4) in the expectation value dominates the integral. This is indeed what 
we see from Figure 3(d), and is consistent with the discussion following (4.6). 
In Figure 3 we also display the effective potential and perturbative ground 
state for the three spacetimes in Figure 1. We see from Figure 2(b) that the 
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bound is not very useful for considering background spacetimes with m 7^ 
and 'large' values of k, as it gives a positive upper bound for unstable solutions. 



5 Discussion 

The existence of perturbatively stable, nakedly singular spacetimes is inter- 
esting, and perhaps surprising if one is a proponent of the cosmic censorship 
hypothesis. We should make some cautionary remarks though, since we have in 
no way demonstrated a violation of this hypothesis. To begin with, one really 
should examine non-spherical perturbations of these spacetimes along the lines 
of the stability proof of the Schwarzschild solution [7,11] and developed further 
in [5]. The existence of non-spherically symmetric unstable modes would indi- 
cate that the existence of any asphericity in a collapsing system would drive 
it away from the nakedly singular configuration, and would not be accessible 
from the study of the spherically symmetric collapse problem. 

However if there continues to be an absence of unstable modes, then we are 
in a somewhat uncomfortable situation. We know from the work of Christo- 
doulou [3] that naked singularity formation in the collapse of massless scalar 
fields is possible, however only with highly tuned initial data. It is likely then 
that massive fields can also form naked singularities, since the behaviour of 
the fields at the singular point is independent of the scalar field potential. If 
the stable branch can be reached from the collapse of suitably regular initial 
data, then we would naively expect that there are choices of initial data that 
end up 'near' the nakedly singular spacetime, and thus the initial data is less 
special in this case. At the very least this is motivation to investigate this 
system further. 

In addition to the solutions considered herein, there are also those with van- 
ishing and negative asymptotic Schwarzschild mass with the same behaviour 
at the centre of symmetry, as well as those with non-vanishing k. These will 
be considered in more detail in the future. 
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A Wyman's Solution 



Here we include some analytic results on exact solutions with V[<f>] = at- 
tributed to Wyman [13]. The equation for the scalar field (2.1c) is integrated to 
find (choosing here k > so that 0o is asymptotically positive and decreasing 
to zero) 

d x [<j>o] = -ke~ 1 2^- x ^/x 2 , (A.l) 

and integrating R 00 = exp(z/ - \ ){dl[v }/2 + (d x [u ]) 2 /A - 9 x [i/ ]9 x [A ]/4 + 
dx[vo]/x) = 0, results in (introducing the integration constant a) 

d x [v ] = ake~^ {uo ~ Xo) /x 2 = -ad x [<i>o], which leads to u = -afa, (A.2) 

where a possible constant of integration has been set equal to zero by requiring 
that {u ,(j) } as x — > oo. Assuming the asymptotic forms of u and 
Ao corresponding to positive and negative asymptotic Schwarzschild masses 
u ~ Tl/ X an( i ~ ±l/x, we find that a = ±l/k. 

Introducing Z := ±k exp(— u /2 + X /2)/(2x) we derive the "transition equa- 
tion" 

xd x [Z] = Z(Z 2 - 2Z/k - 1) = Z(Z - Z+)(Z - Z_), (A.3) 

where Z± := (1 ± VTT¥)/k and we will make use of n± := Z±(Z± — Z T ) = 
±2(a/1 + k 2 ± 1 ± k 2 ) /k 2 . From its definition we require that Z — > as x — > 
from below (above), and from (A.3) we see that Z — > Z± as x — > for the 
positive and negative asymptotic Schwarzschild mass solutions respectively. 
The implicit solution to (A.3) with integration constant chosen so that Z — > 
±k/ (2x) as x — > oo is 

x = ±A(i _ Z /Z + )^+(l - Z/Z.) 1 ^. (A.4) 

In terms of Z we have 

e V0 = ^2(1 + 2Z/fc - Z 2 ), e A ° = (1 + 2Z/A; - Z 2 ), (A.5a) 

h " 'l + 2Z/k- Z 2 )], 



= TA;ln L4^ 

and the large-x behaviour of (A.4) we have 



(A.5b) 



Z~±±.(l±l/x + o(l/x)), (A.6) 
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which is consistent with the required behaviour exp(z/ ) ~ l^fl/x and exp(Ao) ~ 
1 ± 1/x. The small- X expansion gives (defining x k '■= x/y/l + k 2 ) 

Z~Z ± {l-n ± x n k ±), (A.7) 

which leads via (A. 5) to 

e" -^" 2 , e*>~n 2 ± xF, <f>o ~ Tk(n ± - 2) \n(x k ). (A.8a) 
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